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ABSTRACT 

Combination frequencies are observed in the Fourier spectra of pulsating DA and 
DB white dwarfs. They appear at sums and differences of frequencies associated with 
the stehar gravity-modes. Brickhill (1992) proposed that the combination frequencies 
result from mixing of the eigenmode signals as the surface convection zone varying in 
depth when undergoing pulsation. This depth changes cause time-dependent thermal 
impedance, which mix different harmonic frequencies in the light curve. 

Following Brickhill's proposal, we developed analytical expressions to describe 
the amplitudes and phases of these combination frequencies. The parameters that 
appear in these expressions are: the depth of the stellar convection zone when at 
rest, the sensitivity of this depth towards changes in stellar effective temperature, 
the inclination angle of the stellar pulsation axis with respect to the line of sight, 
and lastly, the spherical degrees of the eigenmodes involved in the mixing. Adopting 
reasonable values for these parameters, we apply our expressions to a DA and a DB 
variable white dwarf. We find reasonable agreement between theory and observation, 
though some discrepancies remain unexplained. We show that it is possible to identify 
the spherical degrees of the pulsation modes using the combination frequencies. 

Key words: stars-oscillations;white dwarfs;convection;waves 



1 INTRODUCTION 

Long sequences of almost uninterrupted light curves, ob- 
tained during the Whole Earth Telescope (WET) cam- 
paign on the helium variable white dwarf GD358 (Winget 
et al. 1994, for a more recent analysis, see Vuille et al. 2000), 
disclosed not only the presence of a large number of stellar 
pulsational modes in the Fourier spectra, but also the pres- 
ence of 'combination frequencies', signals that lie at the sum 
or difference frequencies of the stellar eigen-modes (see Fig. 
7 in the first paper). 

Combination frequencies have been observed in other 
pulsating white dwarfs with either hydrogen or helium atmo- 
spheres, e.g., ZZ Psc (aka G29-38, McGra w 1976; Kleinman 
et al. 1998), GD154 ( [Robinson et al. 1978| ), BPM31594 (Mc- 
Graw 1976; O'Donoghue , Warner & Cropper 1992), G117 - 
B15A ( [Kepler et al. 1982| ) and CD 165 ( [Bergeron et al. 1993[ ). 
Indeed, every variable hydrogen white dwarf (class name ZZ 
Ceti, or DAV) that has been observed with sufficiently high 
signal-to-noise ratio exhibits combination frequencies (Bras- 
sard, Fontaine & Wesemael 1995). The same likely holds for 
helium variables (DBV). 

Combination frequencies are thought to result from 
nonlinear mixing of sinusoidal signals that are associated 



with the eigenmodes (named the 'principal modes' in this ar- 
ticle). This conclusion is based on the following arguments: 
combination frequencies are too numerous to be eigenmodes 
themselves (Winget et al. 1994); amplitudes of the combi- 
nation frequencies have been shown to correlate with those 
of their principal modes (for an early review, see McGraw 
1978); combination frequencies tend to have more compli- 
cated fine structure than their principal modes, which can be 
explained naturally by a linear super position of the princ ipal 
modes' rotationally split multiplets (Winget et al. 1994). 

Brickhill (1992) showed that nonlinear mixing arises 



naturally in the context of his theory of convective driv- 
ing (Brickhill 1983, 1990, 1991a, 1991b). He realized that 
the convective turn-over time scale in DA and DB vari- 
able white dwarfs is much shorter than the pulsation pe- 
riod. Thus one can safely assume that the surface convec- 
tive region adjusts instantaneously during pulsation. Brick- 
hill found that under this assumption, the photospheric flux 
variation is delayed and reduced relative to that entering 
the bottom of the convection zone, by an amount depending 
on the depth of the convection zone. Instantaneous adjust- 
ment of the convection zone also implies that the extent of 
the convection zone varies during the pulsation cycle, thus 
leading to variations in the reduction and delay of the flux 
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variation. This distorts the shape of the light curve at the 
photosphere, and brings about the combination frequencies 
in the Fourier power spectrum. Using a numerical analysis, 
Brickhill found that for reasonable amplitudes of the prin- 
cipal modes, he could reproduce the observed amplitudes of 
the combination frequencies. Note that in this theory, the 
combination frequencies reflect distortion of the light curve 
by the nonlinear medium; they are not associated with phys- 
ical displacements and velocities.^ This is indeed confirmed 
observationally (van Kerkwijk, Clemens & Wu 1999). 

In this paper, we use a perturbative analysis to derive 
analytical formulae for the strength and phase of the com- 
bination frequencies. The advantage of this analysis over 
Brickhill's numerical approach is that the dependence on 
stellar properties becomes explicit. We find that two param- 
eters, namely, the depth of the surface convection zone when 
the star is at rest, and the sensitivity of this depth towards 
changes in stellar effective temperature, determine the effi- 
ciency of the mixing process. We also show that two geo- 
metric factors, the spherical degrees of the principal modes, 
and the inclination angle of the stellar pulsation axis, enter 
the analytical expressions. We compare our formulae with 
data on GD358 (a DBV) and G29-38 (a DAV), adopting 
appropriate values for the above stellar parameters. Despite 
imperfect agreement, we show that it is possible to infer the 
spherical degree for the principal modes. In §4.2, we briefiy 
discuss the prospects of explaining the combination frequen- 
cies in other types of variable stars. 



2 PERTURBATION ANALYSIS 
2.1 Origin 

In this section, we demonstrate that the surface convection 
zone can nonlinearly mix sinusoidal fiux variations to pro- 
duce combination frequencies. 

We adopt the following three simplifying assumptions. 
Firstly, we assume neighbouring angular directions do not 
affect each other. Observed pulsations in white dwarfs are 
associated with eigenmodes of low spherical degree = 1 or 
2). Horizontal variations in these modes occur over a scale 
of order the stellar radius. One can safely ignore the effect 
of horizontal heat transfer. This assumption allows us to 
study one angular direction and later generalize the results 
to other directions without modification. 

Secondly, we restrict our analysis to the surface con- 
vective region and assume that the radiative interior caused 
little nonlinearity in the pulsation signals. In the deep ra- 
diative interior where the pulsation is adiabatic, the rela- 
tive importance of the nonlinearity is measured by the local 
pressure perturbation, (Sp/p), which is smaller than a few 
percent for all observed modes. In the radiative region im- 
mediately below the convection zone, the situation is less 
clear. However, it seems reasonable to assume that the non- 
linearity in this region is weak compared to that arising from 
the convective region as the reaction time of the former re- 
gion (local thermal time) is much longer than that in the 



* Fast convection enforces uniform mo vement throughout the 
convective region ( joldrcich fc Wu 1999| ) . There is no distortion 
to the velocity signal. 
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Figure 1. Entropy as a function of pressure in the upper atmo- 
sphere of three hydrogen white dwarf models with similar effective 



temperatures (Wu & Goldreich 1999). The bottom of the convec- 
tion zone (marked with short vertical bars) shifts by one pressure 
scale- height when the surface flux varies by ~ 4%. There is a 
fast convergence in entropy towards the radiative interior. The 
entropy here is dimensionless and is related to the physical one 
by a factor of kg/mp, with kg the Boltzmann constant, and nip 
the proton mass. 



latter (see later). Further study is necessary to confirm or 
to refute this assumption. This assumption allows us to take 
the flux perturbation incident upon the bottom of the con- 
vection zone, {SF/F)b, to be sinusoidal. This also allows us 
to ignore entropy variations in the radiative region (unless 
the material becomes convective). 

Thirdly, we adopt equilibrium models that are adjacent 
in effective temperature to quantify the time-dependent na- 
ture of the convective region. The three models in Figure ^ 
resemble a radial column of gas at flux maximum, passing 
through the rest point, and at flux minimum, respectively. 
This simpliflcation is possible because the thin convection 
zone in a pulsating white dwarf has eddy turn-over times 
much shorter than the pulsation period and can react in- 
stantaneously to pulsation. The entropy of the convection 
zone varies in phase with the surface flux perturbation (this 
is an important feature for driving the gravity-modes, see 
BrickhiU 1983, 1990; Goldreich & Wu 1999, hereafter Paper 
I), the size of the convection zone adjusts accordingly. In this 
analysis, we consider changes in the depth of the convection 
zone caused by entropy variations only. Other effects, e.g., 
convective overshoot and shear turbulence, are not taken 
into account. 

Under these assumptions, we derive the flux pertur- 
bation emerging at the photosphere, {SF/F)pi^, for given 
{5F/F)b. These two flux perturbations are related through 
energy conservation, 

/5F\ ^ /5£\ IdAQ 

VFy'ph \FJt F dt ' ^ ' 
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where AQ is the amount of excess heat (compared to that 
at equihbrium) stored in the convective region, and F is the 
stellar flux at equilibrium. 

Denoting the depth of the convection zone at time t as 
Zb{t), we find 



AQ 



b(t) 



dz pT—{s 



so 



Fn„As + — 
ruv 



dz pT{s — So) 



(2) 



where the second term on the right hand side is small com- 
pared to the first term. Here, all quantities marked with 
subscript are equilibrium quantities, and the thermal time 
constant rt at depth zt is defined as (see Paper I), 

IkB 

F rrir, 



f dzpT. 



(3) 



This time constant is closely related to the conventional 
thermal relaxation time, the latter being defined as rth = 
1/ F{kB /mp) J^'' dz pTcp-, in an isentropic, ionised hydrogen 
plasma, « 2/5 rth- 

The first term on the right-hand side of equation (^ 
quantifies the heat absorbed by the convective region above 
Zbg when its entropy rises uniformly by As. This entropy 
variation is constant throughout most of the region because 
convection carries most of the stellar fiux and because the 
reaction of the convection towards pulsation is fast (Paper I, 
also see Fig. |l]) . The second term in equation (^) represents 
the heat required in expanding or evaporating the convection 
zone. This term is essential for introducing nonlinearity into 
the flux variation. 

Discarding terms higher than second order in As, we 
convert equation (|^) into, 

/SF\ dAs , , , ,dAs 



(-) 



dzbjt) 



, , dAs 



So J 



(4) 



We use s — So = at Zhif) as the entropy perturbation in 
the region z > Zi,(t) is assumed to be zero. Equation (^ is 
similar to equations (40) and (42) in Paper I, except that in 
the present case Th is time-dependent. 

The entropy variation. As, is intimately related to the 
photospheric flux perturbation (eqs. [38]- [39] in Paper I) as. 



As = (B -f C) 



(5) 



where the dimensionless numbers B and C quantify respec- 
tively the response of the photosphere and the superadi- 
abatic region towards gravity-mode pulsation (eqs. [25] & 
[36] in Paper I). In ZZ Ceti stars, B and C are both of order 
8. Again, we take (B + C) to be {B + C){t). 

We define the following dimensionless derivatives, 

91n(B-fC) _ 1 ain(B_+C) 
" 4 91nTeff ' 
91nrb 



/3 



7 



d\uF 
9 In Tf, 1 



c»lnF 4ainrcff' 
and combine equations (^, (j^ and (^ to produce 



(6) 



f)^^X + ..„[l + (2, + ,)X]f . 



(7) 



Ifere, Tc = (-B + C)Ti,, and X is the simplification of 
{5F/ F)p\i. In the temperature range of ZZ Ceti variables, 
our mixing length models yield f} ~ 1.2 and 7 ~ —15 (ob- 
tained using Fig. 1 of Wu & Goldreich [1999]). Note that 
2/9 -f 7 < 0; the thermal content of the convection zone 
increases with decreasing effective temperature. The magni- 
tude of the nonlinear term |(2/3 + 7)^! falls not far below 
unity when X is of order a few percent. Mathematically, 
equation (^) also describes the motion of a pendulum with 
a velocity-dependent mass under the influences of a periodic 
external force and viscosity. 



2.2 Solutions 

We solve for {5F/F)p\^ (= X) using equation (^, first for 
the case when the input signal is a single mode and then for 
the case of two modes. 



2.2.1 Solution for a Single Mode 

At any point on the stellar surface, for a single sinusoidal 
input of the form 



{^^^ = ^» cos (cjit + *0> 



(8) 



the solution for {5F/F)ph can be expanded into 

^ Oi cos {ujit + tpi) + a2i cos {2iJit + tp2i) + (9) 

V ^ / ph 

with the amplitudes and phases obtained from equation (Q) 
as, 



a2i — 



and 



i>2i 



a1 |2/3 + 7l(2^»r,„) 
4 ^1 + (2c^,r,J2 ' 



*i - arctan((Jirco), 
1 



(10) 



= 2il}i + arctan 



(11) 



The expressions for and are identical to equation (66) 
in Paper I; the surface fiux perturbation is reduced in mag- 
nitude and delayed in phase relative to that entering the 
convection zone, as a result of the heat retention (or re- 
lease) of the convection zone during pulsation. Typically in 
pulsating white dwarfs, ranges from a few tenths of a per- 
cent to a few percent. The fact that the harmonic frequency 
of a such a weak mode is actually observable is largely be- 
cause 1 2/3 -f 7 1 3> 1. The amplitude of the second harmonic 
(asi), not considered here, is of order |/37|ai , and should be 
detectable when Oi is large. 

The above solution is best understood in terms of light 
curve distortion. The phase lag and the reduction factor be- 
tween {SF/F)ph and {5F/F)b scale with the thickness of the 
convection zone. Since the convection zone is at its thinnest 
when [SF/F)pYi approaches its maximum, one expects the 
phase lag and the reduction factor to be smaller at the 
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maximum of {5F/F)p\^ than at the minimum. This leads to 
peaked hght curves with sharp ascent and shallow descent. 
Fourier transform of such a light curve exhibit h armonics 
that lead the fundamentals in time, V'2i — 2i/;i > (Brickhill 

199C] 



2.2.2 Solution for Two Modes 

When {SF/F)b is comprised of two sinusoidal signals with 
radian frequencies uii and Uj , 

(^^^ ^ Ai cos {LJit + "iii) + Aj cos {ujt + <iij), (12) 
the following form of solution is adopted, 



V F Jpi, 



ai cos {oJit + tpi) + a2i cos {2ujit + V'2i) 

+aj cos {iUjt + ipj) + a2j cos [2LUjt + il'2j) 

+ai^j cos {{cji ~ ujj)t + ijji-j) 

+ai+j cos {{iOi + ujj)t + ipi+j) + ■•• • (13) 

Subscripts for the different coefficients are chosen to repre- 
sent their corresponding frequencies. We obtain the follow- 
ing generalized expressions for the amplitudes and phases of 
the combination frequencies (including both harmonics and 
mixed combinations), 



Uij Uiaj \2/J + ■-(\{LUi ±LUj)rcg 

1 



ipi±j = {ipi ± ^j) + arctan 



(14) 



(15) 



where Uij counts the number of possible permutations given 
i and j: ntj — 2 if i ^ j, and 1 if otherwise. 



2.3 Angular Integration and Bolometric 
Corrections 

Equations (|l^-([l^) quantify the amplitudes and the phases 
for a combination frequency at every point on the stellar 
surface. In this section, we relate them to the observable 
quantities. 

Defining O and $ to be the spherical coordinates in the 
stellar rotating frame, we adopt the following form for the 
flux perturbation incident upon the convective bottom. 



(^) ^ = ^ (e, $) cos(a;,f 4- m,$ + ) . 



(16) 



Here, p™e""* = Yim, the latter being the spherical har- 
monic function normalized to unity over the sphere. The 'i'i 
appearing in equation (^) is now mi$ -I- ^'^p, while ^'i^ is a 
constant over angle and time. Similarly, the linear part of 
the photospheric flux variation can be written as 



with ai and tpio related to Ai and "^i^ as in equations (|10|) 
An observer with a line-of-sight inclined relative to the 



rotation axis by an angle would detect the following 
bolometric flux variation 



ai5"'(Oo) cos{ujit + ipio), 



(18) 



where the factor g"^ includes effects such as geometric pro- 
jection and limb-darkening. In Table |^, we present the ex- 
pression for g as well as its values when £ < 2. 

The angular dependence of a combination frequency is 
described by the product of the angular dependences of its 
principal modes. This arises because equation ( |l^ is valid 
for every point on the stellar surface. We express the corre- 
sponding disc-integrated flux variations as 

y j ^^ai±jGZ''tj"'HQo)cos{{uJi±ujj)t+tlJig±ja).{19) 

The expression for the G function is presented in the ap- 
pendix, together with some useful values of G. 

Broad-band photometric observations produce ampli- 
tude variations that are related to the bolometric variations 
as (Sf/f) — a\{Sf / f)hoi. For ZZ Ceti stars observed in V- 
band, av ~ 0.4. 

The theoretical expectation value for Rc, the ratio be- 
tween the observed amplitude of a combination and the ob- 
served amplitudes of its principal modes (van Kerkwijk et 
al. 1999), is therefore. 



Rc 



\2f3 + 'y\{i^i±ujj)Tco 



G 



^3 J'^'co) 



(20) 



3 WHAT CAN BE LEARNED? 

In this section, we discuss the potential of extracting infor- 
mation from measurements of the combination frequencies. 
We first describe what information may be available, and 
then compare the observations of two variable white dwarfs 
(a DA and a DB) with our analytical results. 



3.1 Prelude 

A major difficulty of white dwarf asteroseimology lies in our 
inability to securely identify the spherical degree {£) for the 
pulsation modes. In the cases where combination frequen- 
cies are detected, one could use the observed values of Rc 
(eq. 1^) to determine the £ value for the principal modes. 
To illustrate this possibility, we consider the harmonic of 
a m = principal mode. The ratio Gti/lg^gf) (and con- 
sequently the value of Rc) is significantly higher for £ = 2 
than for £ = 1, except when Oo approaches 90° (see Fig. 
^ and eq. [AS]). This arises as the apparent amplitudes of 
higher £ modes suffer stronger cancellation when integrated 
over the stellar disc while the harmonics of these modes do 
not. Note that this is purely a geometric argument and a 
similar method of £ identification would work not only for 



Degeneracy in <& implies that we can take $o = 0. 
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The two dimensionless numbers, /3 and 7, quantify the 
deepening of the convection zone when a white dwarf cools. 
It is therefore related to the width of the white dwarf in- 
stability strip. Let us associate the blue edge of the ZZ Ceti 
instability strip (Toff ^ 12, 000 K) with Tco = 20 s (when the 
lowest order £ — 1 gravity-mode mode satisfies ujtco — 1, 
see Paper I), and the red edge of the strip with Teg — 1300 s 
(when the 1000 s period mode becomes invisible at the sur- 
face, i^Tcg — 10 ^ 1, see Paper I). We find the width of the 
instability strip to be ~ 1000 K when we adopt /3 + 7 ~ — 14 
A larger value of |/3-|-7| would correspond to a nar- 



as in S 2.1 



0„ (deg) 



rower instability strip. These numbers can be obtained from 
combination frequency measurements together with other 
unknown quantities. 

A number of practical difficulties may arise in the actual 
analysis. For instance, different m components of a gravity- 
mode are closely spaced in frequency and may not be re- 
solved by observations of short duration, whereas in obser- 
vations of sufficiently long duration temporal changes in the 
amplitudes of pulsation may occur. In the following sections, 
we apply our results ignoring these difficulties. 



Figure 2. The ratio /{g^ ) as a function of the incli- 

nation angle 0o when £i, Ij take 1 or 2. Notice that at certain 
angles, the combination frequencies may be observable while un- 
favourable geometrical cancellation renders their principal modes 
invisible. 



the sum or difference combinations of two principal modes 
in white dwarfs, but also for other variable stars that ex- 
hibit combination frequencies and where the amplitude of a 
combination frequency satisfies ai±j oc a^aj (as in eq. p^) 
for every point on the stellar surface. 

Inside the pulsational instability strip, the thermal time 
constant of the convection zone (tcq) varies monotonically 
and sensitively with the stellar effective temperature. The 
value of TcQ discloses the relative location of a variable in the 
strip. In addition, we can study convection under the white 
dwarf environment if we can empirically determine the T^s- 
TcQ relation. Time-resolved spectroscop y provides one way to 



measure Tcq (van Kerkwijk et al. 199!:). However, this tech- 



nique requires large telescopes and works only for relatively 
bright white dwarfs. What about using the combination fre- 
quencies? 

The relative phase between a combination frequency 
and its principal modes {tpiQ±jo — {ipio ± ''Pjo)) yields Tcq 
straightforwardly (eq. [|l^). This phase depends on Teg more 
sensitively at low frequency. However, care needs to be taken 
to avoid systematic effects that affect phase measurements 
adversely, such as the presence of small neighbouring peri- 
odicities that are not accounted for. 

Another way to measure Tcq is to take the ratio between 
the amplitudes of the sum and the difference combinations 
from the same pair of principal modes, (eq. ||l^) the geo- 
metric factor cancels when one or both of rrii, mw is 0. 



2^2 G. 



.2 ^"li-^j 
=0 f,- 



(21) 



Note that amplitudes measured in the lower frequency re- 
gion are generally less accurate due to higher noise levels. 



3.2 The DB variable GD358 

For our analysis of the DB varia ble GD358, we use the Whole 
Earth Telescope (WET) data JWinget et al. 1994| ) in which 
different m components of the principal modes are well re- 
solved. We assume that mode amplitudes do not change ap- 
preciably during the entire observation. 



3.2.1 Fine Structure in the Combination Frequencies 

The sum (or difference) combination of two £ = 1 principal 
modes (each split into three m components) can contain up 
to nine fine-splitting comp onents that are different in fre- 
quency (Winget et al. 1994). The relative strengths among 
these components reflect their respective projection in the 
direction of our line-of-sight. In Table |l| we compare the ob- 
served and the theoretically expected values of these relative 
strengths for the sum combination of two principal modes, 
fc = 13 and = 15 (see Fig. 8 of Winget et al. 1994). We 
find overall agreements when the inclination angle Oo falls 
within the range between 40° and 50°. Some mismatches 
exist. They may be caused by the fact that some combina- 
tion components are too close in frequency to be resolved by 
the WET run. In addition, the amplitudes of the principal 
modes vary during the run and this may affect the compar- 
ison. We notice that the WET run lasted 11 days which is 
of order the natural growth ti me for the = 13 and A: = 15 
modes (|Wu fc Goldreich 1999|). 



3.2.2 Combination Frequencies at Large 

Winget et al. (1994, Table 3) determined pulsation power for 
the strongest component in each combination. We compare 
their data with theoretical estimates in Figure ^ For our 
analysis, we assume that all principal modes have spheri- 
cal degree £ = 1 and that the strongest component in each 
combination is produced by the m = components in the 
principal mod es. T he second assumption is tested in Fig- 
ure |. As in ^.2.4 We adopt Go = 45°. We further take 
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principal modes combination frequencies 

fc = 13 fc = 15 V observed amplitudes theoretical 

i (mma) rrij (mma) (mHz) (mma) relative relative 



-1 


6.28 


-1 


9.33 


3033.07 


< 0.5 


< 0.23 


0.23 





5.78 


-1 


9.33 


3038.65 


< 1.0 


< 0.5 


0.21 


-1 


6.28 





19.03 


3039.07 


1.55 


0.72 


0.46 





5.78 





19.03 


3044.65 


2.14 


1.00 


1.00 


+1 


5.46 


-1 


9.33 


3044.76 


< 1.5 


< 0.7 


0.73 


-1 


6.28 


+1 


9.27 


3045.84 


1.79 


0.84 


0.84 


+1 


5.46 





19.03 


3050.76 


1.10 


0.50 


0.40 





5.78 


+1 


9.27 


3051.42 


0.77 


0.36 


0.21 


+1 


5.46 


+1 


9.27 


3057.53 


0.92 


0.43 


0.19 



Table 1. Pulsation ampli tudes for different components in the (k = 13) + (fc = 15) combination. The observed values are extracted from 
Figure 8 of Winget et al. (1994), where the unit mma stands for milli-modulation amplitude and equals to 0.1%. The amplitude of the 
(rrii = —1) + {rrij = —1) component falls below detection limit. The {rrii = 0) + (mj = —1) and (mi = +1) + {rrij = —1) components 
are too close to other components to be resolved. We give rough upper limits to their amplitudes. All relative amplitudes are amplitudes 
scaled by the (m^ = 0) + (mj = 0) component at ^ = 3044.5/iHz. Here, v is the angular frequency and is related to ct) by a; = 2itv. 
The theoretical relative amplitudes are calculated using equation ( po| ) with ©o = 45°. We find that the data are roughly consistent with 
inclination angles between 40° and 50°. 



[2/3 + 7|/qv = 25 and Tcq = 300s to produce theoretical es- 
timates. Such a choice of Teg ensures that aU gravity-modes 
observed in GD358 satisfy the necessary condition for over- 



stabihty, ujtcq > 1 (Goldreich & Wu 199!:). Varying Tcq from 



10 s to 1000 s only changes the estimates by ~ 40%. 

Theoretical estimates based on the above choices of pa- 
rameters largely reproduce the observed values. Figure ^ 
shows that most (~ 90%) combinations that are expected to 
have amplitudes above the ~ 1 mma detection limit are in- 
deed observed. However, a few significant discrepancies are 
present and they merit some discussions. 

The power at = 2660.84^iHz (labelled with 'A' in the 
upper panel of Fig. ^) is attributed to the (fc — 18)-|-(fc = 15) 
combination by Winget et al. (1994). The observed am- 
plitude is roughly twice the theoretical one. Interestingly 
enough, the {k = 16) + (A: = 17) combination lies at 
v = 2659.43/1 Hz and is expected to reach approximately 
the same amplitude as the former combination. This may 
explain the mismatch seen at 'A' if the two combinations 
are not well resolved from each other. 

The combination at 2848.28/iHz ('B') lies ~ 6/iHz away 
from the harmonic of the fc = 15, m = mode and is 
possibly the sum combination between this mode and the 
k = 15, m = —2 mode. As the = 15, m = mode has the 
greatest amplitude in the Fourier spectrum, it is surprising 
that we do not observe its harmonicj^ This problem is not 
unique to GD 358, as our investigation of ZZ Psc (§3.3) finds. 



Full numerical simulations of the convective response may 
provide a solution to this problem (Ising & Koester 1999). 

We can not explain the disagreement seen at 
2946.65/iHz (case 'C') either. 

Combinations at low frequencies are potentially most 



rewarding for estimating Teg . However, the signal-to-noise 
ratio is lower at these frequencies. 



3.3 The DA variable ZZ Psc 

ZZ Psc (aka G29-38) is a relatively cool DA variab le exhibit- 
ing la rge pulsation amplitudes that vary in time ( Kleinman 
1995). We investigate the combination frequencies in this 



star using a set of time-res olved spectroscopy data taken 
with the Keck H telescope (van Kerkwijk et al. 1999). The 
observation lasted five hours and could not resolve the ro- 
tational splitting. However, the signal-to-noise ratio is high 
and phases of pulsation are measured from the data. These 
can be compared with equation (|l5|). 

Among the handful of eigenmodes seen in the spectra, 
one (marked as 'F4') stands out as an ^ = 2 mode (Clemens 
et al. 1999), while all others have £ = 1. This result is sup- 
ported by the fact that 'F4' is associated with a relatively 
large line-of-sight velocity (van Kerkwijk et al. 1999). In this 
section, we show that the combination frequencies also pro- 
vide confirmation for this ^-identification. 

We first assume all observed principal modes are £ — 1, 
m = 0. Taking t^^ = 250 s,^ 12/3 + j\ ^ 10 and Go = 30° ,|^ 
we produce a theoretical amplitude spectrum for the com- 
bination frequencies (Fig. ^), onto which we plot the ob- 
served amplitudes as well as their error bars. The compar- 
ison is satisfactory; all combinations that are estimated to 
lie above the noise level are indeed detected. However, it is 
a surprise to find the three lowest frequency combinations 
to have much higher amplitudes than expected. As signals 
at very low frequencies may suffer from larger noise, this 
discrepancy needs to be confirmed by future observations. 



t Even though fc = 15 is the strongest mode, its harmonic is 
expected to have lower amplitude than, for instance, the combi- 
nation (fc = 15) -I- (fc = 17) due to the factor n^j in equation 
0- 



§ This choice of Tcq ensures that all observed g-modes satisfy 
'^TcQ > 1, and it puts ZZ Psc close to the red edge of the DA 
instability strip. 

^ In this case, we can not constrain |2/3 + 7I and 0o separately. 
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Figure 3. Comparison between observed (triangles) and theoret- 
ical amplitudes (vertical lines) for the combination frequencies in 
GD 358. We calculate all possible two-mode combinations for the 
m = gravity-modes listed in Table 2 of Winget et al. (1994). 
A triangle is solid if the observed frequency of the combination 
lies within 1/iHz of its expected value, where 1/iHz is roughly 
the frequency resolution for this WET run. With this we test the 
assumption that a combination is produced by m = principal 
modes. The top panel expands the central region of the lower 
panel, while the horizontal dashed line in the lower panel is the 
detection limit in this WET observation. Three letters ('A','B' 
and 'C') denote three cases of serious disagreement. 



In Figure y, the combinations noted by 'A' and 'D' are 
respectively the first harmonics of the second largest and the 
largest g- modes. Equation ( |l^ under-predicts the amplitude 
for 'A', and over-predicts it for 'D'. The latter problem, in- 
terestingly enough, appears in GD358 as well (see Fig. |^. 
This may indicate that our perturbation analysis fails at 
large mode amplitudes. The combination at 'C (sum of the 
strongest and the second strongest modes) falling much be- 
low the theoretical estimate jibes with this suggestion. 

Clemens et al. (1999) and van Kerkwijk et al. (1999) ar- 
gued that mode 'F4' has spherical degree £ — 2. If this is 
true, at 6o = 30°, its harmonic would hav e an amplitude 6 
times larger than if it was £ = 1 (see ^3.1 and Fig. |^. This 
could explain the disagreement at 'B'. In the lower panel of 
Figure ^ we show theoretical Rc values for combinations of 
£ = 1 and 2 principal modes. The observed Rc value for 'B' 
is consistent with it being the harmonic of an = 2 mode. 
This is a piece of independent evidence supporting the £ = 2 
identification for mode 'F4'. At 0o = 30°, Rc values for com- 
binations that involve 'F4' and another mode do not differ 
appreciably from other combinations (see Fig. ^). 

In the upper panel of Figure |^, we study the phase 
difference between a combination and its principal modes. 



Figure 4. Amplitudes of combination frequencies in ZZ Psc. The 
observed values are plotted as solid triangles with error bars. The 
vertical lines are theoretical amplitudes of all possible combina- 
tions between observed gravity-modes. They are crowned with a 
open circle if they lie within 50/^Hz in frequency from the ob- 
served combinations. This is roughly the frequency resolution for 
this set of data. The noise level in the data is indicated by the 
horizontal dashed line. 



From equation ([l5|), one expects most of these phase differ- 
ences to be positive and to lie close to 0. Indeed, the observed 
phase differences cluster closely around 0. However, in detail 
the fit is not good. This may be related to the unaccounted- 
for periodicities which influence the phase determinations 
(van Kerkwijk et al. 1999) 

Employing data from WE T and other observations (see 
Kleinman et al. 1998), Vuille (2000b) analyzed the relative 
phases of combination frequencies in ZZ Psc. He found 
tpi±j — {tpi ± i/'j) to be small but predominantly positive. 
This reflects the fact that pulsation light-curves have sharper 
ascent than descent. Vuille concluded that these relative 
phases remain fairly constant over the years, a finding con- 
sistent with equation (p^). 



4 CONCLUSIONS 

4.1 DA & DB Variables 

Our analysis leads to two key formulae (eq. and p^) 
that describes the strength and the phase of the combination 
frequencies relative to their principal modes. 

A few stellar parameters enter these formulae. They are 
the thermal constant of the stellar convection zone at equi- 
librium (tcq), the rate of deepening of the convection zone 
with cooling of the star (quantified by \2f3 + 7I), and the 
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Figure 5. For ZZ Psc, we compare observed values of relative 
phase and amplitude ratio between combination frequencies and 
their principal modes, plotted as triangles, with theoretical expec- 
tations. Relative phase (eq. Jist) does not depend on the spherical 
harmonic degrees of the two principal modes, while amplitude ra- 
tio (Re in eq. |2o| ) does. We plot Rc for 1=1 and 2 principal 
modes. The combination with the largest Rc value is the har- 
monic of mode 'F4', while the other two circled combinations are 
mixing results of 'F4' with other principal modes. 



Combination frequencies are produced by the surface 
convection zone in a pulsating white dwarf. Photosphere in 
these stars is not thermally capable of distorting the light 
curve. We therefore expect equations and (|l5|) to hold 
for all wavelengths. 



4.2 Other Types of Variables 

Combination frequencies have also been reported in two 
PG1159 variables (PG1707-h427, Fontaine et al. 1991; 
HS2324+3944, Silvotti et al. 1999). Presumably, these hot 
white dwarfs do not have surface convection zones. What 
could be distorting the light curves? 

Could a radiative, partially ionising layer produce the 
distortions? Such a layer is believed to exist in the upper 
atmosphere of PG1159 variables and is believed to be re- 
sponsible for driving the observed pulsations. It is similar 
to the surface convection zones in DA and DB variables in 
that it retains heat when warmer, and releases heat when 
cooler. However, unlike in the case of the convection zones, 
the amount of heat retained (or released) by the partial ion- 
ising region can not be significantly modulated throughout 
the pulsation cycle by the presence of other pulsation modes. 
This is because the reaction time of the ionising region is 
roughly the local thermal relaxation time, which is of the 
same order as periods of overstable modes. Thus, we can 
not explain the combination frequencies in PC 1159 stars. 

We note that in other types of small amplitude pul- 
sators, e.g., 5-Scuti stars, sdB variables, and 7-Doradus 
stars, combination frequencies have also been reported. We 
conjecture that a thin surface convection zone is present in 
these variables and is capable of exciting pulsation modes, 
as well as distorting the light curves. 



inclination angle between the observer's line of sight and 
the stellar pulsation axis (6o)- It is becoming possible to 
use the combination frequencies to constrain these stellar 
parameters. We find that for both GD358 and ZZ Psc, the 
observed amplitude spectra can be roughly reproduced us- 
ing reeisonable choices of the above parameters. The same 
choices can also explain the values for the dimensionless 
numbers (a, 6 and c) Brickhill (1992) summarized from his 
numerical study. The I and m values of the eigenmodes also 
enter into these formulae. This presents the potential of de- 
termining the I values of the principal modes using the com- 
bination frequencies. An 1 = 2 mode is expected to have a 
stronger harmonic than an £ = 1 mode, and this is indeed 
observed in ZZ Psc. 

When analyzing observed combination frequencies, we 
have ignored amplitude variability during a long observing 
run, or in the case of a short run, have assumed that modes 
are axisymmetric with respect to the pulsation axis (m = 0). 
The failure of these assumptions may account for some of 
the discrepancies between observation and theory and may 
prevent us from accurately determining stellar parameters. 
More suitable data sets might yield more conclusive infor- 
mation. 

We find that theory over-predicts the amplitude in the 
harmonic of the strongest pulsation mode in the two stars 
we considered. We suspect that it results from the stronger 
nonlinearity associated with the largest mode. 



The author would like to acknowledge the many ben- 
eficial comments and suggestions by Drs. Peter Goldreich, 
Joerg Ising, Marten van Kerkwijk, Scot Kleinman and Fran- 
cois Vuille. 
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^ \ \m\ 
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- - 0.28 

1 - 0.245 sin 00 0.346 cos 00 

2 0.126 sin2 0o 0.126 sin(20o) 0.05[1 + 3 cos{20o)] 

Table 2. Values of gV^ as functions of 0o, the inclination angle between the observer's line-of-sight and the stellar rotation axis. 
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Table 3. Values of G™'^™^ /(g™' (?™^ ) as functions of 0o. Values of G™' lidg.^ dg^') a^rs obtained by reversing the sign of rrij 
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APPENDIX A: ANGULAR INTEGRATION 

Let (0, $) be the spherical coordinate system defined by 
the stellar rotation axis, and (0, (/>) be that defined by the 
observer's line-of-sight. Let 0o be the angle between this 
line-of-sight and the rotation axis. The two coordinate sys- 
tems are related by 

cosB = — sin Oo sin 6 cos </> -I- cos Oo cos ^, 
sin cos $ = cos Oo sin cos (ji + sin Oo cos 0, 
sinOsin"I> = sin S sin;/). (Al) 

To obtain the photometric pulsation amplitude, we in- 
tegrate the photospheric flux variation over the visible disc, 



1 

2^ 



MM)MdM(^) (e,$,t), (A2) 



where [i — cos 9, and h{^) is the limb-darkening function 
normalized by h{fj,)fj,dfi — 1. For our exercise, we adopt 
the Eddington limb-darkening law of /i(/i) = 1-1- 3/2fi (see, 
e.g., Shu 1991), as is appropriate for grey-atmosphere. 
It is convenient to define 



We list values of g in Table Q for < 2. 

The angular dependence of a combination frequency is 
described by the product of the angular dependences of its 
principal modes. Its photometric amplitude is related to its 
photospheric amplitude by the following function, 

'^Z'tj'"^ = 27r '^'^ J ^^^^^'^^^ 

xP™'(e, $)P,"^' (6, $) cos((m, ± m,)$). (A4) 

It is easiest to evaluate the G function by reducing the above 
product of the Legendre functions into a linear sum of such 
functions. We present values of the ratio Gi.\, ^ / (g^"' 3f . ^ ) 

in Table ^ for £i — £j — 1. Values of /"^ / (g™'(?™^ ) can 
be trivially obtained from the sa me t able by changing the 
sign of rrij. For our discussion in §3.1, we also need 



~tO±0 

,-rO±0 
^2 2 

9°29°2 



1.97 + 0.82 cos-^ Qo 
Scos^ Oo — 1 

8.39 + 2.51 cos(2eo) + 0.22 cos(4eo) 
(3cos2 Go - 1)2 



(A5) 



d(t> / h(n)ndnRe[Yr{e,^)]. 



(A3) 
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